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Abstract 

These lecture notes address an audience of physicists or mathemati- 
cians who have been exposed to a first course in quantum mechanics. We 
start with a brief discussion of the general "system-bath" paradigm of 
quantum dissipative systems, analyze in some detail the simplest example 
of "pure dephasing" of a two-level system, and review the basic concept of 
the density matrix. We then treat the general dissipative time-evolution, 
introducing completely positive maps, their relation to entanglement the- 
ory, and their Kraus decomposition. Restricting ourselves to Markovian 
evolution, we discuss the Lindblad form of master equations. The notes 
conclude with an overview of topics of current interest that go beyond 
Lindblad Markov master equations. 

These notes were prepared for lectures delivered by F. Marquardt in Octo- 
ber 2007 at the Langeoog workshop of the SFB/TR 12, "Symmetries and 
Universality in Mesoscopic Systems" 

1 Introduction 

The following general situation is of interest in many fields of quantum physics, 
ranging from quantum optics to condensed matter physics: A single quantum 
system interacts with a large reservoir, alternatively called "bath" or "environ- 
ment" . Whenever the system is driven out of equilibrium by external perturba- 
tions, this coupling makes the system relax back to equilibrium. 

The system itself might be a single atom, the spin degree of freedom of an 
electron or a nucleus, a single particle moving through some potential landscape 
(such as in a man-made interferometer in a semiconductor), or even a many- 
particle system. In most of these cases the system has only a few degrees of 



freedom, and sometimes it even has a finite-dimensional Hilbert space (such 
as for the spin). In contrast, the bath needs to have infinitely many degrees 
of freedom in order to generate truly irreversible dynamics. It might be the 
electromagnetic field into which the atom can radiate its energy, or the crystal 
lattice that gets distorted by an electron moving along. 

In a classical setting, the bath introduces 
dissipation (friction), as the system's energy 
can irreversibly be transferred to the bath. 
By necessity, this also means that the system 
will experience a fluctuating force. The bal- 
ance of those two effects makes the system 
settle into thermal equilibrium at a temper- 
ature set by the bath. In quantum dynam- 
ics, there is yet another feature beyond these Figure 1: A quantum system in- 
two: The system can display coherent effects, teracting with the infinitely many 
i.e. interference phenomena in space or time. degrees of freedom of an environ- 
These will be destroyed gradually by the cou- ment. 

pling to the environment, an effect known as "decoherence" or "dephasing". 

In the following lecture notes, we will start with a simple example that in- 
troduces the basic physics and requires only knowledge of elementary quantum 
mechanics, applied to the two-level system. It will be used in subsequent sec- 
tions to illustrate the concepts. We then review the description of "mixed" 
(incoherent) quantum states by way of the density matrix. Afterwards, we dis- 
cuss the most general framework for describing dissipative quantum dynamics: 
completely positive maps. In the next step, we restrict ourselves to the impor- 
tant and simple class of Markov dynamics and introduce the Lindblad master 
equation, which is a workhorse of open systems dynamics in many branches of 
physics. Finally, we outline a few topics of modern research that go beyond this 
elementary tool. 




2 Example: Pure dephasing of a two-level sys- 
tem 

2.1 Effect of a classical stochastic process 

We consider a system described by a two-dimensional Hilbert space, Ti = C 2 . 
This is a two-level system, called "qubit" in some modern applications. Since 
the spin 1/2 is the most important physical realization, we denote the two basis 
states as "spin up" ||) and "spin down" |J.). The vector describing the system's 
state contains the two complex probability amplitudes for these basis states: 




In the absence of a magnetic field (i.e. in the absence of an energy difference 
between the two spin states), ipyi(t) would be time-independent. We now look 
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at the simplest possible model for dephasing of this two-level system: The energy 
difference between both basis states fluctuates in the course of time. For the 
moment, we will not even keep the bath as a large quantum system in its own 
right. Rather, we assume that the fluctuations it produces can be treated as a 
classical stochastic process. This becomes a good approximation for the limit 
of high temperatures. The time-dependent Hamiltonian we thus want to study 
is 

*(.)-(*>!!)• (2 > 

Here e(t) is the stochastic process, which has the meaning of a fluctuating mag- 
netic field along the z-direction if we think of the case of a spin. The Schrodinger 
equation ih^ \ip(t)) = H \ip{t)) can be solved directly. The components of the 
solution are 



V-iW = ^(0) (3) 
Mt) = e*<*>VT(0), (4) 

where the fluctuating phase is the integral over e(t): 

< p {t) = -\j\{t>)dt>. (5) 

I 1 2 

First of all, we observe that the probabilities V'T/I W do n °t evolve, and only 
the relative phase between both states is affected by the noise. This is the 
reason for speaking of "pure" dephasing in this example. Next, we turn to look 
at observables that are sensitive to the relative phase between the two basis 
states, e.g. the operator o x that is proportional to the spin's component in the 
x-direction: 

Its quantum-mechanical expectation value, for a fixed realization of the process 
e(i), is given by 

m) \a x \m) = ^I(0)V T (0)e^W +c.c. (7) 

[Remarks on notation: We employ standard "bra-ket" notation, (tp \A\ 4>) = 
(ip\A(j)} = ^A(f>, where {ip \ip) = 1. The complex conjugated version of the first 
term on the right-hand-side is abbreviated as c.c] 

In reality it will be necessary to repeat the experiment many times in order 
to obtain an estimate for the expectation value, by averaging over the outcomes 
of a x measurements (which yield either +1 or —1 in each individual run). In 
general, each run will correspond to a different realization of t(t) as well. There- 
fore, the result obtained by this averaging procedure actually also involves taking 
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Figure 2: (a) Pure dephasing can be viewed as arising from different, fluctuating pre- 
cession frequencies for an ensemble of spins; (b) Frequency fluctuations as a stochastic 
process; (c) Resulting phase fluctuations; (d) Visibility 



the classical (stochastic) expectation value of the random phase factor e llf '. We 
will denote this average by (•) in order to distinguish it from the quantum 
mechanical expectation value. Thus, the experimentalist will observe 

{W)\°*W))) V = W^T(O) (e^'^ + c.c. (8) 
Apparently, the factor 

v{t) = (e 1 ^ (9) 

describes the modification of the "interference term" ■0*(O)V't(O) due to the 
noise. The same factor would occur if we were to evaluate, say, the expecta- 
tion value of a y . We will here call this factor "visibility", as it has to do with 
the contrast or visibility of interference experiments that rely on the coherence 
between the two basis states. In particular, the overall magnitude of the inter- 
ference term is suppressed by \v(t)\ < 1. This will usually decrease in the course 
of time, though that is not guaranteed and will depend on the properties of the 
process e(i). 

If e is a Gaussian stochastic process (of zero mean), and therefore tp(t) is a 
Gaussian random variable with (ip(t)} — 0, we can express the visibility explic- 
itly in terms of the variance of the phase: 

v(t) = e-^ v{t)2 ^. (10) 

2.2 Several ways to obtain the same mixed state 

If we are only concerned with the system and the expectation values of any of its 
observables (i.e. operators acting on the system alone), there are indeed many 
more situations that will lead to the same predictions. Let us enumerate them: 

(a) The influence of a classical stochastic process, with the subsequent aver- 
age, has been treated above. 

(b) Interaction with a bath in a random initial state: Suppose we are con- 
sidering a large Hilbert space, consisting of system and bath: Hsb = Us 

The product states of the form \tps) ® \xb) G 7~Lsb (where \ips) S an d 
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\Xb) € Hb) will be denoted as \ips) \xb) for brevity, where it is understood 
that the first part refers to the system, and the second to the bath. 

In order to describe the interaction between system and bath we might cither 
write down the full Hamiltonian, or else (more conveniently for our purposes 
here) state the action of the full unitary time-evolution operator Usb^) that 
maps an initial state in Hsb onto a final state at time t. Given an orthonormal 
basis of states in Hb, we postulate that the time-evolution simply results 
in a phase factor that depends on the state, in the form: 

UsB(t) |T) \Xj) = e^ (t >|T)|x,) (11) 
UsB(t)\i)\xs) = |i)bG>. (12) 

This defines U$b- Now suppose furthermore that the initial state of the bath 
is picked out of the \Xj) at random, each of them occuring with a probability 

Wj (where J2j w j = !)• This means that the full initial state is V'sbW^ ~ 
\ip(Q)) \Xj) with probability Wj. Such a situation actually arises if the bath is 
at a finite temperature, where different energy eigenstates occur with classical 
probabilities fixed by the Boltzmann distribution. 

We now want to calculate the observed average value of a x at time t. To 
this end, we must first evaluate the quantum-mechanical expectation value 
of the product operator (<r x ) s (g> (1) B , with respect to the state V'ssW) = 

UsB(t) iPsb{0)^, an( i then average according to the weights Wj. We write the 
operator as a x for short, it being understood that a x acts only on the system 
and leaves the bath's state untouched. Straightforward calculation using the 
rules pT| and ( fl2| shows that we obtain exactly the same result as before: 

SS {^skt) \°M { p B {t)) = wm(o) ( eMt) ) v + c - c -< ( 13 ) 

3 

if we identify (e 1 ^^} with Wje lV} [the extension to a continuous probabil- 
ity distribution is obvious]. Moreover, the results are the same for any system 
operator whose expectation value we care to evaluate. In this sense, the inco- 
herent ( "mixed" ; see below for definition) system states produced according to 
(a) and (b) are indistinguishable. 

(c) "Quantum randomness" : Instead of postulating classical randomness in 
the choice of the initial bath state, we could as well have postulated that it is a 
superposition of the basis states, 

|x(0)>=£V%IXi>> (14) 

3 

which is normalized since . (y/uij) 2 = 1. Under the action of UsB(t), the 
initial state 

\i> SB (Q)) = W> T (o) IT) + ^|(0) II)) (8 |x(o)) (15) 
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will evolve into 



\^ SB (t)) = v T (o) IT) |xtC*)> + V^i(o) II) \xi(t)) ■ (16) 

Here we have defined the two bath states IxtW) an d |X|(*)) that have evolved 
out of |x(0)) under the influence of the system being in state ||) or \l), respec- 
tively: 

3 

\xi(t)) = EV^Ix;> = lx(o)>. (18) 

3 

Again, when calculating the expectation value of a x , we obtain the same result 
as above. However, in addition we realize that the visibility v(t) can now be 
expressed as the overlap of the two bath states |xt/J.(^)) : 

( e ' VW ) = EX e ^ (t) = <*i(*) IXt(*)> • (19) 

3 

This is a result that is much more general than the present example would 
suggest. The two states are sometimes called "pointer states", referring to the 
quantum theory of measurement. As it were, the environment can be thought 
of as measuring the state of the system, which means that |x(0)) evolves into 
either |xt(t)) or |xi(i)), depending on whether the system had been in the state 
|t) or |J,). This is reminiscent of the pointer of a measuring device pointing 
either way depending on the signal it picks up. When the two pointer states 
become orthogonal, this implies that perfect knowledge about the system could 
be inferred (in principle) from the state of the bath. As a consequence, inter- 
ference effects in the system itself are completely destroyed, analogous to the 
gedanken experiment of the "Heisenberg microscope" [3] • 

In the most general framework, one can set up a path-integral analysis, 
where the bath's state would evolve according to the particular trajectory of the 



system. The overlap ( 19 1 between two states having evolved under the influence 
of two different system trajectories is then called the "Fcynman- Vernon influence 
functional" [7] . This forms the starting point for path-integral evaluations of the 
dissipative dynamics of some systems, such as the damped quantum harmonic 
oscillator or tunneling decay under the influence of dissipation. For further 
details see the book of Weiss ITBl. 



3 The density matrix 

The examples of the previous subsection all yield the same results for the ex- 
pectation values of any system operator. It is therefore useful to introduce a 
description of the system's state that makes this equivalence explicit and which 
is general enough to deduce from it any arbitrary expectation value. This is 
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achieved by the density matrix, whose properties we review in this section. It 
forms the main object of study in the field of open quantum systems, where 
the goal will usually be to calculate the time-evolution of the density matrix 
of a system subject to a fluctuating environment. While originally the word 
"state of a quantum system" referred only to "wave functions" , i.e. vectors in 
a Hilbert space, it is now extended to include the "mixed states" that have to 
be described by a density matrix. Most readers will know this material and can 



skip directly to subsection 3.4 which returns to the example of pure dephasing. 



3.1 Classical uncertainty 

There are different ways in which the density matrix concept may arise. First, 
let us assume there is some classical uncertainty about the system's state, i.e. 
the state G Tis occurs with probability Wj. As explained above, when 
trying to predict the experimentally observed average value of an observable 
A, we thus should first take quantum-mechanical expectation values and then 
perform an additional, "classical" average: 

(A) = |>W (20) 

3 

We can obtain the same result by introducing the density matrix as a weighted 
sum over projectors \ipj) (4>j\, 

p=Y,wA^)m, (2i) 

3 

and taking the trace over pA: 

(A)=tr(pA). (22) 

Here p S B(H), where B(TL) denotes the set of bounded operators on the Hilbert 
space TL. 

In general, this step (from the set of \ipj) to p) involves a compression of 
information. Regardless of the number of states \ifij) we started out with in 
the beginning (which need not be orthogonal!), p itself is an operator on Tis, 
represented by a dmiHg x dim7i5 matrix in the case of a finite-dimensional 
Hilbert space. 

3.2 General properties 



The following properties follow directly from the definition (21 ) 



trp = 1 (23) 
pt = p (24) 
P > (25) 
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For a "pure" state \ip), the density matrix is a projector onto that state: 
p = \tp) (tp\. As a consequence, p 2 = p and trp 2 = trp = 1. In general, the 
"purity" trp 2 obeys trp 2 < 1. All states that are not pure are called "mixed", 
and their density matrix has a rank larger than one: 



trp 2 < 1 4^ pis mixed. (26) 

When making a measurement that is able to distinguish the particular state 
\<p) from other, orthogonal states, we can phrase this by saying that we measure 
the operator \<f>) (<f>\, which will yield the value 1 only when we indeed find \<p). 
Thus, the probability of observing \<f>) is the expectation value of its projector: 



P* = tr(p|0)M) = ^|p|0. (27) 

In the special case where p describes the pure state \ip), this correctly reduces 
to the well-known postulate of quantum mechanics, P ^ = \(cf> \ip)\ 2 - 

After diagonalizing p, its eigenvalues pj can therefore be interpreted as the 
probabilities of finding the respective eigenvectors in a measurement that distin- 
guishes between those eigenvectors. The properties p > and trp = 1 are then 
seen to correspond to the simple fact that these probabilities are non-negative 
and normalized. 



3.3 Reduced density matrix 

We now turn to a different setting in which the density matrix arises: Consider 
the world made up of a system and a bath, as explained in the introduction. For 
generality, we will describe the overall state of the world by a density matrix p SB 
(understood in the sense explained above). Suppose we are now only interested 



in evaluating expectation values of system operators As- Adapting (22 1, we 
have 

(A s ) = (A s <8> 1 B ) = tr SB (p SB (A s <8> 1 B )) . (28) 

Here we have indicated that the trace is taken in the full Hilbert space TLsb- 
However, since the operator acts only on the system alone, we can break this 
trace into two steps: 

(A s ) = tr s (tr B (p SB )A s ) = tr s {p s A s ) . (29) 

Here we have introduced the reduced density matrix of the system, by taking 
the partial trace over the bath degrees of freedom: 

p s =tr B (p SB ) (30) 

We can make this explicit by choosing a product basis in TLsb- Then the matrix 
elements of p s are given by 

Pi'i — ^ Pfi^.iXi.i) > (31) 
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where refers to the basis state \ipi) \xj), with £ Hs and \xj) G Hb 



It is straightforward to confirm that the general properties (23), (24 1, and 
(25), continue to hold for p s . 



3.4 Example: Pure dephasing 

For the example described in section [2] the density matrix is given by 

m= ( IV> T (o)l 2 ^(0)^(0)«(i) \ 

This is independent of the specific model, i.e. it holds regardless of whether we 
think of a classical noise process or of the interaction with a quantum- mechanical 
environment (where p would be the reduced system density matrix). Only 
the off-diagonal element p n (t) = (| \p(t)\ J.) = ^ T (0)^(0)«(t) is affected by 
dephasing. The decay of \v(t) \ will reduce the purity of the state [we abbreviate 

p=hMci)l 2 ]: 

tr(p 2 (t)) = l-2p(l-p)(l~\v(t)\ 2 ) (33) 

When starting from an equal-weight superposition, p — 1/2, this tends to 1/2 
if v — > 0. In that case, one ends up in the fully mixed state: p = 1/2. 

Let us have a look at the simplest possible case for the decay of v, postulating 
that it decays exponentially in time: 

v(t) = (e l{p ^ =e- r * t , (34) 

where is called the "dephasing rate". Such a decay will occur when the 
underlying fluctuations of e are of the "white noise" type, and therefore (p(t) 
undergoes Brownian motion {ip is proportional to a Wiener process). Then 
(ip 2 (t))^ oc t, which leads to exponential decay according to (exp[iip(t)]) v = 
exp[— {ip(t) 2 ^ /2] = exp[— T v t]. For this particular case, we can write down a 
simple first-order differential equation for the elements of p(t) : 

f>u(t) = -r^PnW ptt = pu = °- ( 35 ) 

This is the simplest example of a "Markov master equation" . The general 
structure of such an equation, which we will discuss later, is of the form 

p = hp. (36) 

However, it is not permissible to postulate an arbitrary operator L on phe- 
nomenological grounds, because that might turn out to violate the basic prop- 
erties of p. This can be seen clearly in the following example: Suppose we want 
to describe the spontaneous decay of the excited state of an atom by emission of 
a photon. The first reasonable (and indeed correct) ansatz that comes to mind 
would be to postulate an exponential decay of the probability p^ of finding the 
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atom in the excited state: p^ = —Tp^, and consequently p\n = —p^ = +Tp^, 
to conserve probability. Suppose we were to assume that the off-diagonal ele- 
ments do not change (pn = p^ = 0), which seems to be the simplest possible 
assumption. In the long-time limit t — > oo, this would lead to 

„(t- oo) = ( m ° (0) Pn } 0) ) (wrong). (37) 

This is not a positive scmidcfinite matrix if pn(0) ^ 0. In other words, such an 
ansatz would violate a basic requirement, the positivity of probabilities. The rest 
of these lecture notes is concerned with describing the correct genera il structure 
of the time evolution of density matrices, which makes sure that such pitfalls 
are avoided. 

4 Completely positive maps and Kraus opera- 
tors 

4.1 Complete positivity 

Let us consider the linear map that takes the initial density matrix p(0) to its 
value at time t, 

$ : B(H) -» B(H) 

p(0) ~ p(t). (38) 

Later on we will consider the map $ for different times, but for now we suppress 
the corresponding subscript for brevity. 

Let us list the properties which we will require of <i>, which follow from the 
properties of p that we want to be respected by the time-evolution: 

tr<f>(p) = tip (norm — conserving) (39) 
<&(p') = $(p) (preserves hermiticity) (40) 
p>0 => $(p)>0 ("$ is a positive map") (41) 

It would seem that these three features are all that is needed to have a permis- 
sible time-evolution. Surprisingly, that is not the case. The third requirement, 
of <& being a positive map, has to be replaced by the stronger condition of <& 
being "completely positive". 

Definition- We call a linear map <& : B(H) — > B(H) 11 completely positive" 
(CP) iff the following holds: Given any other Hilbert space H', we consider the 
product space H = H (£> H' and construct a linear map : B(H) — > B(H) 
that derives from $ by acting like <!> onto the original Hilbert space Ti and 
not affecting the space H', that is: <&(/? ® p') = <f>(p) ® p' for p e B(H) and 
p' G B(H'). Then is a positive map. 
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4.2 Example of a positive but not completely positive map 
(relation to entanglement theory) 

At first sight, it is surprising that positivity of $ itself is not enough to guarantee 
the positivity of the enlarged map We will now look at the simplest example 
of a map that is positive but not completely positive. This example is of interest 
in the theory of entanglement. 

Claim- <& : p i — > p* is positive but not CP. 

It is obvious that $ is positive, as transposition does not change the eigen- 
values of a matrix. Now let us consider the enlarged map operating on a 
product basis in TL = TL ® TL' . It induces a "partial transposition" of a density 
matrix p £ B(TL): 

Mp)](i',j')(.hi) = hij'WJ) ' ( 42 ) 
where we observe the interchange of indices i and i' referring to Ti. In order 
to prove that $ is not a positive map, it is enough to find one example of a 
state p > (for one suitably chosen Hilbert space TL') for which the partial 
transposition fails to make Q(p) > 0. It is clear that any product state, of the 
type p = p® p' , will not be sufficient for this purpose, since p > and therefore 
®(p) > as well. More generally, due to the linearity of the same holds for 
any so-called "separable" state which is a mixture of product states: 

P = £ w jP ij) ® P' U) ■ ( 43 ) 

3 

Here p^',p'^ are valid density matrices in B(7i) and B(Ti'), respectively, and 
Wj > 0, Wj = 1. Any separable state will yield a positive semi-definite par- 
tial transpose: Q(p) > 0. This is the important PPT ("positive partial trans- 
pose") criterion, a necessary condition for separability of a state p, discovered 
by Asher Peres in 1996 [15] . 

In order to find an example without a positive partial tranpose, we thus have 
to consider non-separable, i.e. so-called " entangled" states. The simplest case 
is having both Ti and TL' two-dimensional Hilbert spaces, and considering the 



pure state p 



with 



^) = ^(ITT> + IU». (44) 



The partial transposition acts like $(|TT) (III) = lit) (Til (note the interchange 
in the first position, referring to TL), and likewise on other combinations occuring 
in the projector p. As a result, we find 



Hp) 



^ 1 ^ 

10 

10 

\ 1 J 



(45) 
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where the matrix has been written down with respect to the basis 

ITT) > ITT) i ITT) j ITT)- It has eigenvalues +1/2 (three-fold degenerate) and —1/2. 

A a consequence, <& is not positive, and $ is not a CP map. 

As a side-note, we mention that PPT is even a necessary and sufficient 
criterion for separability if dim7Y = dim7Y' = 2 as in our example [and this 
even holds when the dimensions are 2 and 3, respectively]. One has to go to 
higher-dimensional Hilbert spaces in order to find states that are not separable 
(i.e. entangled) but still have a positive partial transpose. 

Physically, the property of complete positivity ensures that one ends up with 
a permissible state even if the dissipative system had been entangled initially 
with another system. As a consequence, there is no analogon to the concept of 
complete positivity for classical dissipative systems, since classical physics does 
not know entangled states. 



4.3 Kraus decompositions of a CP map 

It turns out that all CP maps that fulfill the other properties mentioned above 
can be decomposed in a simple way. We state the theorem, due to Kraus [TTj . 
without proof: 

Theorem- Provided we are given a map $ that fulfills the properties d39| 



and (40 1, as well as complete positivity, there exists a set of "Kraus operators" 
Kj : B{TL) — > B{TL) that are normalized in the sense Y2j KjKj = 1 and that 
can be used to represent <&(/o): 

®{p)=Y J K jP K}. (46) 

3 

The converse also holds (i.e. the three properties follow from the existence of 
such a decomposition). In the case of an infinite-dimensional Hilbert space, the 
set of values for the index j can be countably infinite. 

We note that this decomposition is not at all unique. For example, multi- 
plying Kj by a phase factor e l9j changes nothing. More generally, a unitary 
matrix with matrix elements Uj/j can be used to convert to a new set of Kraus 
operators that represent the same map: K'y = Yl j KjUjji [where the cardinality 
of the two sets is assumed to be the same, adding Kraus operators equal to zero 
if need be] . It can be shown |16j that any two equivalent Kraus decompositions 
are connected in this way. We also note that for a finite-dimensional Hilbert 
space any map can be represented using at most (dim7i) 2 Kraus operators. 



4.4 Examples 

We now list a few examples for the applications of Kraus decompositions: 

• Purely unitary time-evolution: \tp(t)) = U(t) | -0(0)) leads to 

p(t) = U(t)p(0)U' (t), and p(t) — <J>(p(0)) can therefore be represented by 
the single Kraus operator K — U(t). 
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Random unitary evolution: If we pick some Uj at random with prob- 
ability Wj, as in the example with pure dephasing by classical noise, 
then the set Kj = ^/w]Uj(t) will yield the correct p(t) = <f>(p(0)) = 
J2j w jUj(t)p(0)Uj (t) , and normalization follows from J2j w j = 1- 

However, in the example of pure dephasing of a two-level system, we can 
also choose a more economical decomposition. For example, for v(t) e R, 
just two Kraus operators suffice: K\ = 1^/(1 + v)/2 and 
K 2 = <r zy /(l-v)/2. 

The a z Kraus operator in the previous example can be thought of as 
describing a "phase flip" (changing the sign of the off-diagonal clement 
of the density matrix). In the same sense, u x would describe a "bit flip" 
(turning |f) into |J.) and vice versa), and a y a combination of the two. 

In quantum information processing, two-level systems are viewed as "quan- 
tum bits" (qubits) . After sending such a qubit through a communication 
channel (where it can be subject to technical noise, or even interact with, 
and get entangled with, some bath), its state will have been changed by a 
map <& that is characteristic for this quantum channel. Thus a quantum 
channel can be described by giving a set of Kraus operators. 

As indicated by the two- level example from above, for a finite-dimensional 
Hilbcrt space a finite number of Kraus operators will suffice (even if the 
underlying physical description of the noise involved an infinite number of 
possible different time evolutions Uj ) . This is very important in the theory 
of quantum error correction. It means that, contrary to first appearances, 
quantum computers are not as bad as classical analog computers when it 
comes to error correction. 

Kraus operators can be used to describe measurements: After an ideal 
von-Neumann measurement that distinguishes between the orthogonal 
states \4>j), the system ends up in one of those states with probability 
Pj = tr(p {4>j\)- However, from the point of view of someone who is 
not told the measurement outcomes, the state after the measurement is 
described by the density matrix p' = J2jPj \4>j) (4>jV an d the map from 
p to p' can be described by the Kraus operators Kj = \<j>j) {<j>j\. The 
more general case of measurements that reveal only partial information 
(POVM: positive operator valued measurements) is described by arbitrary 
Kraus operators that are not necessarily projectors, and the probability 
of finding a particular value is then tr (pK^K^j . 

For systems with high- (or infinite-) dimensional Hilbcrt spaces (like those 
with a continuous degree of freedom, e.g. a harmonic oscillator), the Kraus 
decomposition, though still possible in principle, becomes less useful in 
practice due to the large number of Kraus operators. 
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4.5 Construction of Kraus operators for system-bath in- 
teraction 



We now give an explicit construction of the Kraus operators for the impor- 
tant example where the time-evolution of the reduced density matrix is due 
to the interaction with an environment. Suppose that at time t = the en- 
vironment is in the initial state |x(0)), uncorrelated with the arbitrary initial 
system state p s (0). This means for the total (system+bath) density matrix: 
p SB (0) — p s (0) (g |x(0)) (x(O)l- I n addition, we allow for an arbitrary unitary 
time-evolution Usb acting on the product Hubert space TLsb, taking the full 
initial state to the final state at time t. Therefore, the system's reduced density 
matrix at time t is: 



p(t) = tv B (u SBPs ®\x(o))(x(o)\ui 



= T,^\ U sb\x(0))ps(x(0) 



SB 



Xj 



In the second line we have introduced the sum over a basis \xj) of the bath 
Hilbcrt space TLb to perform the trace over the bath. The resulting matrix 
elements of the full time-evolution Usb with respect to the bath states define 
the operators Kj , as indicated. These operators act on the system Hilbert space 
TLs- They are normalized, since 

E X K = E (*(0) y sB \xi) (Xj \Usb\x(0)) = 1. (48) 
3 i 

As a consequence, we can identify them as the Kraus operators needed to de- 
scribe the map of p(0) = p s (0) to p(t). If the bath initially were in a mixed 
state, we would find p(t) to be the weighted average of expressions of this type, 
and consequently the full set of Kraus operators would consist of the individual 
sets, multiplied by the factors where Wj are the weights. In this way, 

the time-evolution due to interaction with an environment (starting from an 
uncorrelated state) can always be expressed using Kraus operators, i.e. as a CP 
map. The idea of the general construction of Kraus decompositions, presented 
in [IB], makes use of this fact. 



5 Markov master equations of Lindblad form 

In the previous section, we have been dealing with the completely arbitrary 
time-evolution of a density matrix. Let us now specialize to evolutions of the 
Markov type, i.e. where the density matrix follows an equation of the form 
p = Lp, with L representing a linear map. This is called Markovian since the 
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evolution of the state only depends on the present state itself, not on its history. 
Such an equation is the direct quantum analogue of the evolution equation for 
the probability density for the case of a classical stochastic process of Markov 
type. 

5.1 Lindblad's theorem 

It is clear that the Markov property, when combined with the restrictions dis- 
cussed in the previous section, will lead to a particular form of L. This problem 
has first been considered in its full generality by Lindblad [13] (1976). In writing 
down the assumptions for his proof, which we will list below, Lindblad considers 
the map in the Heisenberg picture, where the time-evolution is applied to the 
observable A instead of the density matrix p: 

tr[p(t)A] = tr[* t (p(Q))A] = tr[p(Q)$f (A)], (49) 

for any choice of initial density matrix p(0). 

Then a "completely positive dynamical semigroup" is defined as follows: Let 
A be a W^-algebraQ and <I> t be a family of completely-positive maps of A into 
itself, depending on the real- valued time parameter t S M.q . In addition, we 
require the following properties: 



3f(l) = 1 (50) 

*? •*? = $f +t (VM>0) (51) 

lim||$f-l|| = (52) 

<I>^ is normal (ultraweakly continuous) (53) 



t->0 



The first line guarantees conservation of the trace of the density matrix. The 
second line is the central semi-group property, which will guarantee Markov 
dynamics and represents a strong restriction on the allowed physical situations. 
It tends to work as a good approximation in situations where the correlation 
time of the fluctuations characterizing the environment is short in comparison 
with other time scales, such as those set by decay rates and oscillation periods. 
Under these conditions, Lindblad proved the following: 

The action of can be expressed in the form of a Markov master equation, 
where p(t) = <I>t(/?(0)) fulfills the differential equation 



p(t) = Lp(t), (54) 

A W*-algebra (also called a von-Neumann algebra), is a *-algebra of bounded operators 
on a Hilbert space that is closed in the weak * topology and contains the identity operator. 
This essentially means it is an algebra of operators A, which can be added and multiplied 
with a scalar as usual (thus forming a vector space), have the usual rules of non-commutative 
multiplication between operators, the possibility of taking the adjoint (that is what the * 
refers to) with all the features you would expect, and an operator norm that is compatible 
with the adjoint operation: = ||j4|| 2 . See books on functional analysis, such as |17| . 



15 



with the Liouvillian operator of "Lindblad form" : 

Lp = -~ [H, P ] + ^2(RjpR] - \r]RjP - \pR]R 3 ) ■ (55) 

3 

Here the first part, involving some hermitian Hamiltonian H, describes the stan- 
dard unitary time-evolution of p (possibly with renormalized matrix elements 
of H due to the presence of the environment, i.e. H ^ H$, where Hq would be 
the intrinsic Hamiltonian) . The dissipative dynamics is generated by the second 
term, with the relaxation (or Lindblad) operators Rj that do not have to fulfill 
any special constraint (unlike the normalized Kraus operators, from which they 
derive). We note that, just as for the Kraus operators, the choice of Rj is not 
unique, and even the separation into a unitary and a dissipative part is not 
unique either. 



5.2 Obtaining the Lindblad form from the Kraus decom- 
position 

We now give the main ideas of the derivation, building on the Kraus decompo- 
sition, without attempting mathematical rigour. 

The density matrix at a small time 8t deviates from p(t — 0) only to first 
order in St, due to continuity and the Markov structure. In addition, it can be 
written in the form of a Kraus decomposition: 

p(6t) = p(0) + O(St) = ]T K jP (0)K}, (56) 

3 

where Y]j KjKj — 1. In order to satisfy this structure, we need one Kraus 
operator close to unity, which we will write as 

K = l + (-^H + A)6t+ ... (57) 
n 

The second term has been decomposed into "real" and "imaginary" parts, with 
two hermitean operators — H and = A. All the other Kraus operators 



Kj (with j > 1) must be of order 0(\/5~i), to obtain the desired expression (56 1 
We can now obtain A from the normalization requirement: 



which yields 



1 = E K \ K i = 1 + mA + E K s K i . ( 58 ) 

3>0 j>l 



A =-^X K ^ (59) 

3>l 



that is an expression of O(St ). Inserting this result first into (57 1 and then into 



the decomposition (|56|, we obtain the Lindblad structure (551, by identifying 
the relaxation operators as 
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5.3 Examples 

We now list a few examples for the relaxation operators occuring for different 
dissipative processes. In the present lecture notes, we will not address the 
techniques used for obtaining microscopic derivations of these operators and 
the rates occuring in there, which depend on the specific type of environment 
and the assumed coupling. In general, however, the relaxation operators are of 
a simple form, describing the operator that induces the dissipative transition, 
multiplied with the square root of the corresponding rate. 

• Pure dephasing of a two-level system, with exponential decay of the visi - 



/ r 

bility at a rate T v : R = \ -^-a z yields the Markov master equation (35) 



Exponential relaxation from the excited state to the ground state, at a 
rate T: R = \/Ta_, where ct_ = (f|. This yields the proper decay of 
the probabilities, p^ = —Tp^ and p^ = —pii = +Tp^, but it also gives 
non-trivial information on the decay of the off-diagonal elements: 

i r 
Pn = ~^( £ t ~ e i)Pn ~ 2^ti ■ ( 61 ) 

The first term results from the unitary evolution (i.e. from H = ej |T) (T| + 
e l ID (41); whereas the second term describes dephasing at a rate r® ff = 
r/2 that is exactly half the total decay rate. If pure dephasing is present 
in addition, one just needs to introduce a second relaxation operator, as 
given above, and the rates would add: r^ otal = Y v + r/2. In the context 
of nuclear magnetic resonance or qubit physics, the fact that r^ otal > 
r/2 is often expressed as an inequality for the corresponding time-scales 
(inverses of the rates): T 2 < 2Ti, where T 2 = (r^ 1 )" 1 and T x = Y~ x . 
We have encountered this inequality here ultimately as a consequence of 
the requirement that the time-evolution preserve positivity of p [which 
any decay rate smaller than r/2 for the off-diagonal elements would not 
ensure] . 

If the bath is at finite temperature, the two-level system might also absorb 
a quantum of energy, i.e. become thermally excited. This is described by 
R = v/r\^cr + , where er + = al_ = |f) (J.|. 

Taking into account all the three processes discussed up to now leads to 
the so-called "Bloch equations" first invented to describe the dissipation 
and decoherence of systems such as atoms or spins. 

Damping of a harmonic oscillator (due to a coupling to the bath that is 
linear in the coordinate of the oscillator) is described by R — vTa. Here 
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r is the damping rate (which can be observed in the linear response of the 
system), and a is the annihilation operator that reduces the occupation 
number of the oscillator by one: a \n) = yjn \n — 1). [a = (x + i^j)/(2xo) 
with x = yjh/ (2muo); x and p being position and momentum operators 
for the oscillator] There are many important applications: For example, 
a single standing wave mode of the electromagnetic field inside an optical 
or microwave cavity is described as a harmonic oscillator. That oscilla- 
tor is damped because the photons can leak out of the cavity through 
the semi-transparent mirrors of the cavity, and the individual decay pro- 
cess corresponds to the destruction of a single photon inside the cavity: 
tl i ► tl 1 . Another example of current interest arc nanomcchanical sys- 
tems (small beams on the micrometer scale) which are harmonic oscillators 
to a good approximation, damped because of their connection to the me- 
chanical structure onto which they are attached and into which they can 
radiate phonons. 

We have not described how to obtain the rates and the operators, and under 
which conditions one may expect the physical dynamics to be well approximated 
by a Markov master equation. While a thorough discussion of these points would 
go beyond the scope of the present notes, and it is hard to list simple, generally 
valid conditions, the following rulc-of-thumb may be offered: most systems in 
which this approximation works have decay rates that are both much smaller 
than the typical transition frequencies of the system, and also smaller than the 
inverse correlation time of the fluctuating force that the environment imposes 
on the system. This condition is usually achieved in the limit of a weak coupling 
between system and bath, as the decay and dephasing rates become arbitrarily 
small in that limit. 

6 Beyond Lindblad equations 

Although in practice Lindblad Markov master equations are used in the majority 
of applications in the various subfields of physics, current research in quantum 
dissipative systems focuses on the interesting effects that arise outside of this 
framework. Here we just list a few of those physical situations and features, to 
give the reader a sense of what lies beyond Lindblad dynamics. 

• In some areas it is experimentally feasible to measure the quanta which are 
emitted into the environment by the system, thereby learning more about 
the system's state. For example, the photons having been emitted from 
an atom or leaking out of an optical cavity can be registered by photo- 
detectors. The dynamics of the reduced density matrix, conditioned on 
the observed detector results, may then often be described by a modified 
master equation that depends on these results. Such an approach some- 
times is referred to as "quantum jump trajectories simulations" , because 
it was initially developed to describe quantum jumps in individual atoms 
that have been observed by light scattering. 
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Figure 3: Possible features of the visibility in the example of pure dephasing (section 
[1J|, due to non-Markovian dynamics (a) and due to non-Gaussian noise (b). 

• Non-Markovian dynamics is generated when the environment's fluctua- 
tions display long correlation times. A generic consequence of taking into 
account these effects are short-time deviations from exponential decay. In 
addition, under some circumstances the coherence of the system can even 
be revived at later times, after having decayed initially. Non-Markovian 
dynamics is generated when the environment's fluctuations display long 
correlation times. A generic consequence of taking into account these 
effects are short-time deviations from exponential decay. In addition, un- 
der some circumstances the coherence of the system can even be revived 
at later times, after having decayed initially. Non-Markovian dynamics 
is generated when the environment's fluctuations display long correlation 
times. A generic consequence of taking into account these effects are 
short-time deviations from exponential decay. In addition, under some 
circumstances the coherence of the system can even be revived at later 
times, after having decayed initially. 

• A favorite microscopic model for the environment is a bath of harmonic 
oscillators, which is, for example, an exact representation of the electro- 
magnetic field or the phonons in a crystal lattice. When such a bath 
is coupled to a single particle (bilinearly in the coordinates of parti- 
cle and bath oscillators), one speaks of the "Caldeira-Leggett model" 
[3J HI [TH]. This allows to answer questions that go beyond Markov dy- 
namics. For example, when such a bath gives rise to standard diffusive 
motion (<V 2 (i)) oc t) at high temperatures, the "Quantum Brownian mo- 
tion" resulting at zero temperature is sub-diffusive, with the distance to 
the origin obeying (r 2 (i)} oc ln(i). One can use the same kind of model to 
study the suppression of quantum tunneling due to the influence of a dis- 
sipative environment. These nonperturbative studies are usually carried 
out in a path-integral framework. 

• As the coupling to the environment is increased, there can be qualita- 
tive changes in behaviour at some critical coupling strength. The best 
known example occurs in a model where a particle that tunnels between 
two states is coupled to a bath of harmonic oscillators. This is commonly 
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referred to as the "spin-boson model" , where the "spin" refers to the two- 
level system and "boson" refers to the bath of oscillators. For a certain 
distribution of bath oscillator frequencies ("Ohmic bath"), increasing the 
coupling strength beyond some point makes the system undergo a quan- 
tum phase transition towards a symmetry-broken phase where the particle 
remains trapped in one of the two states for all times. Ergodicity is broken 
due to the strong dissipation. 

• Even if the actual physical environment is not a bath of oscillators, it 
often may be treated as such to a very good degree of approximation as 
long as the coupling is weak. The fluctuating force generated by a bath 
of oscillators is Gaussian-distributed. However, for strong coupling, this 
approach will fail, and the dissipative dynamics can become qualitatively 
different due to the non-Gaussian fluctuations acting on the system. An 
example of current interest are the current fluctuations generated by the 
passage of discrete, single electrons through nanostructures. 

7 Further reading 

Mathematical treatments of quantum dissipative dynamics may be found in the 
books by Davies |6J (published in 1976, still without reference to Lindblad, but 
developing the same concepts) and Kraus [TT] , who emphasizes the connections 
with the theory of measurement. The treatise of von Neumann on the mathe- 
matical foundations of quantum mechanics |10j forms the basis. 

A thorough introduction to the density matrix and its uses, at an elementary 
level for the physicist, can be found in the book by Blum [1], which also presents 
a microscopic derivation of the rates and relaxation operators appearing in mas- 
ter equations. Dissipative quantum systems in general (often with emphasis on 
quantum optics) are described in more detail in the books by Carmichael [5], 
Gardiner and Zoller jB] , and Breuer and Petruccione |Hj . The book by Weiss [IB] 
emphasizes those aspects that go beyond the Lindblad master equations, such 
as the Feynman- Vernon influence functional formalism, some exact solutions, 
and a very detailed discussion of the spin-boson model. In that context, the 
classic review by Leggett and co-workers on the spin-boson model is also highly 
recommended [T2"] . 

Kraus operators in the context of quantum information processing are de- 
scribed in the book by Nielsen and Chuang [14 , as well as in the quantum 
information lecture notes by Preskill ([T3], chapter 3), who also describes the 
Lindblad master equation and gives a nice proof of the Kraus representation 
theorem. 
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